Laser-driven atom moving in a multimode cavity: strong enhancement of 

cavity-cooling efficiency 
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Cavity-mediated cooling of the center-of-mass motion of a transversally, coherently pumped atom 
along the axis of a high-Q cavity is studied. The internal dynamics of the atomic dipole strongly 
coupled to the cavity field is treated by a non-perturbative quantum mechanical model, while the 
effect of the cavity on the external motion is described classically in terms of the analytically obtained 
linear friction and diffusion coefficients. Efficient cavity- induced damping is found which leads 
to steady-state temperatures well-below the Doppler limit. We reveal a mathematical symmetry 
between the results here and for a similar system where, instead of the atom, the cavity field is 
pumped. The cooling process is strongly enhanced in a degenerate multimode cavity. Both the 
temperature and the number of scattered photons during the characteristic cooling time exhibits 
a significant reduction with increasing number of modes involved in the dynamics. The residual 
number of spontaneous emissions in a cooling time for large mode degeneracy can reach and even 
drop below the limit of a single photon. 

PACS numbers: 32.80.Pj, 42.50.Vk, 42.50.Lc 
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I. INTRODUCTION 

The mechanical effect of light on atoms in high-Q 
cavities is currently being the subject of intensive the- 
oretical and experimental research. The large inter- 
est that this field attracts is due to the complexity of 
the strongly-coupled dynamics of a moving atom and 
a few-photon field, as compared to the relatively well- 
understood atomic motion in external laser fields. Vari- 
ous effective cooling and trapping schemes have recently 
been found, which can drive atoms to low, sub-Doppler 
temperatures on a fast time scale. The use of dynami- 
cal fields in optical cavities is also a promising approach 
to extend the power of ordinary laser cooling methods 
to a large variety of atomic species. There are cooling 
mechanisms independent of the internal electronic struc- 
ture of the particle, which could enable us to achieve the 
principal goal of cooling molecules. 

In the simplest cavity-cooling schemes jj], ^, [| [|, |j| 
the cavity is used as a passive element to taylor the 
atomic spontaneous emission rate by altering the spec- 
tral mode density. For example, in the pioneering work 
by Mossberg et al. Q , an atom moving in a resonant laser 
standing-wave field can undergo a Sisyphus-type cooling 
in coloured vacuum, e.g. the one created inside a cavity. 
The spontaneous emission is then favoured, by selection 
in frequency space, from the top of a potential hill into 
the bottom of a potential well. Similarly, Vuletic et al. 
H pointed out that a two-photon Doppler-effect can oc- 
cur provided the atom tends to scatter inelastically an 
incoming laser photon into a cavity mode of higher reso- 
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nance frequency. 

Dynamic cavity-cooling effects have been first pre- 
dicted in Refs. || In these schemes the free atom 
moves in a weakly driven cavity field. The damping re- 
lies on a strongly coupled atom-field system, in which the 
field exhibits non-adiabatic dynamics due to the finite 
mode bandwidth k. Let us emphasize that the cooling 
mechanism is quite sophisticated, requiring at least two 
elements: (i) strong atom-field coupling manifested by 
the dressed state spectrum of the atom-cavity system; 
(ii) a cavity bandwidth in a specific narrow parameter 
range for the non-conservative dynamics. Both condi- 
tions are present in current setups J8|, ^ and indications 
of the predicted mechanical effects have been observed in 
experiments. 

In this paper we will consider nonadiabatic effects in 
a configuration similar to the ones in Refs. (I], |[ |[ |), 
that is, a moving atom is driven by a classical laser field 
from the side, and cavity photons are created only via 
the atom, by scattering from the laser field. We pro- 
ceed in two respects. First, we move towards the strong 
coupling regime where the non-conservative dynamical 
effects, analogous to the ones explored in Refs. || 0, ap- 
pear and amount to new cooling regimes in the parameter 
space. Accordingly, we will present a non-perturbative 
quantum treatment of the internal dynamics of the cou- 
pled atom-cavity system. Our model is exact in the low- 
saturation regime that happens either for large atomic 
detuning or for weak external driving intensity. The dif- 
fusion coefficient due to dipole heating and the linear 
friction coefficient will be analytically calculated, which 
allows for the statistical description of the center-of-mass 
(CM) motion characteristics. We obtain the somewhat 
surprising result that pumping the atomic or pumping 
the field component of the coupled system are not the 
same even for low saturation. Although the similarity 
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appears in the form of a quite intuitive symmetry be- 
tween the results we derive and the ones presented in 
Refs. H [?J for the cavity-driven system, important phys- 
ical differences follow from the two different situations. 
Such consequences can be found, for example, in the tem- 
perature dependence on the system parameters. 

As a second step forward, we will study multimode ef- 
fects of the cavity-induced forces. Multimode fields have 
been conjectured to give rise to a significant improvement 
in terms of steady-state temperature when many cavity 
modes are driven simultaneously |l0| , and also in the case 
of the two-photon Doppler effect HT In the physical con- 
figuration we consider in this paper, the only component 
of the system being externally driven is the atom. It is a 
very natural and, in practice, a simple extension to put 
a multimode cavity, such as a confocal resonator, around 
the driven particle. We ask what happens if there are 
several modes interacting with the atomic dipole. Here 
we will concentrate on scaling laws rather than on the 
geometrical aspects of the multimode cavity field, this 
latter is being important, e. g., to reveal atomic trajec- 
tories jl2|]. We consider degenerate modes with closely 
uniform mode functions. This situation proves to be an 
easy way to effectively enhance the dipole interaction 
strength and reach previously unexplored regimes of the 
dynamics. One of the main questions is if the number of 
spontaneously scattered photons during the cooling time 
can be restricted below the single photon level, which is 
the necessary requirement for effectively cooling particles 
without closed two-level pumping cycle. 

The paper is organized as follows. In Sect. 2 we present 
the model of the system, and introduce the definition of 
the mechanical forces and the diffusion that we use to de- 
scribe the semiclassical CM motion. In Sect. 3 the set of 
linearized Heisenberg-Langevin equations is solved which 
leads to the analytical expressions for the linear friction 
and the diffusion coefficients. These results are analyzed 
first for a single-mode cavity in Sect. 4, with the aim of 
identifying the cooling mechanisms relying on a passive 
cavity and those originating from the strongly- coupled 
atom-cavity dynamics. Then, in Sect. 5., we study how 
the physically relevant quantities, such as the steady- 
state temperature, the cooling time, and the number of 
spontaneously scattered photons scale with the number 
of degenerate modes of a multimode cavity. Finally, we 
conclude in Sect. 6. 



II. THE MODEL 

The system is composed of M degenerate modes of a 
high-finesse cavity interacting with a moving atom. The 
atom is driven from the side by a coherent laser field (see 
Fig. Q). We restrict the study of the atomic motion to the 
longitudinal dimension along the cavity axis. The driving 
field can be assumed uniform in the relevant range of the 
atomic motion. This is equivalent to neglecting the slow 
variation of the Gaussian envelope which is broad com- 
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FIG. 1: Schematic representation of the system composed of 
a cavtity and a moving atom that is driven from the side. 



pared to the wavelength. The pumping strength is thus 
given by a single parameter rj. The pumping frequency 
ujp is detuned from the atomic transition frequency with 
A a = LUp—ujA and from the common resonance frequency 
of the modes with Ac = uj p — ujc- The dipole-cavity field 
interaction is described by the Hamiltonian 

M 

H = —hAc a\a.i — TiAa^ cr — ihrj(a — a^) 

i=l 

M 

~ iK^imizW a i ~ 9t(z)a\cj] , (1) 
i=l 

which is written in a frame rotating at the pumping fre- 
quency uip. The atomic dipole operator is denoted by cr, 
and a, is the bosonic field-amplitude operator associated 
with the ith m ode. The p osition-dependent dipole cou- 
pling is <?.; = \JTiloc /2eVi dfi(z), where d is the atomic 
dipole moment, Vi is the effective volume of the mode i, 
and fi(z) is the respective mode function which is nor- 
malized to have the maximum equal to 1. 

In addition to this coherent dynamics, the system is 
subject to dissipation processes, i.e. to spontaneous emis- 
sion with rate 27 into radiation modes other than the cav- 
ity ones, and to cavity photon loss with a uniform rate of 
2k for all the degenerate modes. The resulting open sys- 
tem can be described by Heisenberg-Langevin equations. 
We consider the low-saturation regime where the atomic 
population inversion operator a z = (a^a — aa')/2 can be 
replaced by the c-number —1/2. With this approxima- 
tion one gets a linear set of equations 

M 

cr = (iA A - j)<t - 22gi(z)ai + r\ + £ A (2a) 
i=i 

d, = {iAc - k)<2j + g*(z)a + Ci , (2b) 

for all the modes i = 1..M. The Langevin noise operators 
are defined by the second-order correlations 

(£(iiK f (*2)) = 2 7 *(t 1 -i a ) (3a) 
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all the other correlations vanish. Note that a simple 
transformation, a — ► a — 77/ (iAjy — 7) shifts the driv- 
ing from the atomic dipole a to the mode amplitudes 
Oj, thus one could think that driving the field and driv- 
ing the atomic component leads to equivalent dynamics. 
However, this is not true for a moving atom when z is 
considered also as a variable. 

In the Heisenberg-Langevin equations of the field am- 
plitudes and the atomic dipole, the atomic position ap- 
pears as a parameter. Nevertheless, under the mechani- 
cal effect of the radiation field on the atomic center-of- 
mass motion, the atomic position is also a variable. In 
the present paper we use a consistent semiclassical model 
that relies on a classical treatment of the CM motion [jll| . 
The external atomic CM dynamics is assumed to decou- 
ple from the internal one given by the Eqs. (|^), and it is 
governed by a Langevin-type equation 



x = p/m 

p = f P + Pp + £, p , 



(4a) 
(4b) 



where f p is a conservative force, f3 is a friction coefficient 
and £ p is a classical stochastic term giving rise to dif- 
fusion. Statistical properties of the CM motion can be 
derived from this Langevin-equation. This semiclassical 
approach is consistent as long as the temperature is far 
above the recoil temperature T roc = h 2 k 2 /2m, where k is 
the wave number of the modes involved, m is the mass 
of the atom. 

The aim of the present paper consists in calculating 
the coefficients in Eq. (^b). They all can be derived from 
the force operator, which is given by the gradient of the 
interaction term in the Hamiltonian W\\ 



M 



F = inJ2[(d z 9i)^a i -(d z g*)a\a} 



(5) 



First, the mean of the force operator provides directly f p , 
f P = (F) ■ (6) 

Then, the noise term £ p originates from the fluctuations 
of the force operator. On evaluating the two-time corre- 
lation function, 



(F(h)F(t 2 )) - (F(t 1 ))(F(t 2 ))=V p S(h 



(J) 



the momentum diffusion coefficient can be identified with 
the coefficient of the function S(t\ — £2)- Note that here 
we use a different approach to the diffusion coefficient as 
compared to earlier work §, 0, 0, |l|. It has the ad- 
ditional advantage that the motional diffusion is linked 
directly to the quantum noise accompanying the dissipa- 
tion processes of the system, i.e., the dipole fluctuations 



due to spontaneous emission (£4 that scales with 7) and 
the cavity photon loss that scales with k). It allows 
to reveal how the basic noise sources of the internal dy- 
namics yield diffusion of the external CM motion. 

Finally, for the friction coefficient f3, the dependence of 
the internal variables dj and a on the velocity has to be 
taken into account. Indeed, the non-conservative friction 
force arises from the time-delayed reaction of the internal 
variables to the displacement of the atom. They require a 
period of kT 1 or 7 -1 (in the coupled system typically the 
shorter time, though this depends also on the coupling 
strength) to reach the stationary state adapted to the 
actual position of the atom. If the atom moves much 
less than the wavelength during this time, i.e. in the low- 
velocity limit kv -C 7, k, the velocity-dependence of the 
variables can be taken into account in a consistent way. 
The time derivative is to be replaced by d/dt — > d/dt + 
vd/dx and, simultaneously, the variables are expanded 
into the series 
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The resulting dynamical equations can be solved system- 
atically in different orders of the velocity v. The first 
order components obey the set of linear equations 



(iA A - 7 )a« - J2 9i(z)aP = d z a^ ( 8a ) 

i=l 

(iA c -n)at ) +g*(z)a^ = d z af , (8b) 

The coefficient (3 of the linear friction force can be defined 

as 



A I 



[3/m = <F«) = zh^idzg*)^ } af ] + a^> af) - H.c. 

(9) 



III. SOLVING THE LINEAR SYSTEM 

In this section we perform the calculation of the linear 
friction and the diffusion coefficients by directly solving 
the equations describing the internal dynamics, i.e. Eqs. 
(||) and (||). Let us introduce the Fourier transform of 
the operators, 



0(t) 



1 



/2tt 



dn e - int 



o(n) 



(10) 



and use the rule d/dt — > — iVt. The velocity-independent 
part of the system variables yields 
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t7<°>(fi) 



i 

D(0) 
1 

D(n) 



N 



k=l 



{k - iA c - m)(V2^s(Q)r, + £i(n)) - 5] 



W) ^ ( n)-f;#^(n) 



K — jAp 



iAc — 



fc=l 



K — ll\C — 1 



(11a) 
(lib) 



where the determinant of the linear system reads 



JV 



D{Q) = (k- iA c ~ ifl)(7 - iA A \9k{* 



fe=i 



(12) 



The system variables are obtained as linear combinations 
of the pumping term r\ and noise operators. Note that, 
as a consequence of the coupled dynamics, both the field- 
amplitude operators et^ and the dipole operator a incor- 
porate the fluctuations £4 and associated with the 
spontaneous emission and from the cavity photon loss, 
respectively. 

In any normal-ordered expression of the operators aj, 
cr\ a, and a, all the noise operators £ are found on the 
right side, while all the adjoint operators occur on the 
left side of the expression. Therefore, when evaluating 
the mean value of a normal-ordered product, the noise 
terms have no contribution. It is useful to introduce the 
c-number variables which arise from the coherent driving 
terms of the exact expression (^lj). They correspond to 
the stationary solution of a semiclassical model of the 
same system, and read 



8<® = V 



k — iAc 
D 



I D 



(13a) 
(13b) 



where, and hereafter in the paper, D — D(0). In normal- 
ordered products, the operators can be replaced by these 
simple "semiclassical" solutions. For example, the atomic 
excitation is obtained as 

(at ( °V°)) = fl<°> V°> = 7/V + A 2 C )/\D\ 2 , (14) 

which has to be well below one in order to be consistent 
with the low-saturation assumption. The steady-state 
photon number in the mode i is 

= (^ 0) ^)\ 9l (z)\ 2 /( K 2 + A 2 c ). (15) 



The second expression exhibits how the photon number 
is related to the atomic excitation. As this latter is nec- 
essarily small, the photon number has to be below g/ ' k. 
A. The friction coefficient 

The friction coefficient is defined, in Eq. (^J), by a nor- 
mally ordered expression. Hence, it is enough to take into 
account the coherent driving terms, i.e., the semiclassical 
solution (]l^) of the variables. Accordingly, the compo- 
nents, first-order in velocity, have to be determined from 
the Eqs. (||) by using only the c-number part of the so- 
lutions. This simplifying fact is in accordance with the 
physical intuition. The non-adiabatic dynamics of the in- 
ternal variables, being the origin of the damping, does not 
substantially depend on the noise. On the other hand, 
the non-adiabaticity must be reflected in the dynamics of 
the c-number semiclassical variables, since it includes the 
pumping and damping processes of the coupled internal 
system. 

The first-order components of the variables can be eas- 
ily obtained. They are 



^) = ^[( K - i A c ) 2 -G}(d z G) + ^T 
( K -iA c +j-iA A )g*(d z G) 



(16a) 



„(i) 
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V 



D 2 n-iA 



where 



M 

fc=i 

M 



(D(d z g*)-g*T) , (16b) 



r = X>(0,<£), 



(17a) 
(17b) 



fc=i 



Note that (d z G) = T + T*. Replacing these expressions 
into the definition (JoJ) , one gets the friction coefficient. 
Here we present only the solution for a standing-wave 
cavity, where the mode functions are real. It is 
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{d z Gf 



-Ah 



kA a + 2A C [k + 7 + k 



« 7 - A C A A + G/2 
k 2 + A% 



if 



\D 



(d z G) 2 A c (kA a + 7 Ac)[k 2 A a + 1 2 A c + (A A + A C )(A A A C - G)] 



(18) 



kA( 



D\ 2 n 2 + A 2 

1 ° k 



r 



The generalization for running- wave modes in a ring cav- 
ity is straightforward. 



B. Diffusion coefficient 



general validity regardless of the relationship of the pa- 
rameters involved. The only condition is the weak atomic 
excitation that can always be met by a proper adjustment 
of the pumping strength 77. 



The diffusion stems from the fluctuation of the force 
operator, i.e., that of the dipole interaction term of the 
Hamiltonian (|]J). When calculating the diffusion coeffi- 
cient from its definition (]?]) , the product of the force op- 
erators contains non-normally ordered terms. Hence, the 
noise terms in the solutions (|ll|) have non- vanishing con- 
tributions, and the noise correlations given in @ have to 
be used in the calculation. One can separate terms orig- 
inating from the spontaneous emission noise (I; a) and 
from the cavity loss noise (£i). We omit the details of 
the lengthy calculation here apart frommentioning one 
non-trivial step in the derivation. Although the sources 
£4 and £j (i = 1...M) are supposed to represent white 
noise, the resulting noise associated with the force op- 
erator becomes coloured, that is, it has a non-uniform 
spectrum. Accordingly, in addition to the Dirac-(5 cor- 
relation assumed in the definition (Q), one obtains other 
terms proportional to the derivatives of the Dirac-<5 to 
all order. We neglect these terms and consider the coef- 
ficient of the Dirac-(5 to describe the diffusion process. 

The diffusion coefficient for real mode functions reads 



P dip = 2fi 2 -^ U^(W 



(d z G) 2 A ( 



kA a + 7 A C 
\D\ 2 



(19) 



In addition to the dipole heating, the noise induced by the 
random recoil accompanying the spontaneous emission 
has to be taken into account. The recoil contributes to 
the total diffusion by 

V roc = 2h 2 k 2 A u 2 j^(K 2 + A 2 c ) 1 , (20) 

where k A — lo A /c and u 2 , characteristic of the atomic 
transition, is the mean of the recoil momentum projected 
on the cavity axis. 

To conclude this section let us emphasize that the re- 
sults for the friction and the diffusion coefficients have a 



IV. COOLING REGIMES 

In this section we evaluate the previously calculated 
expressions for the diffusion and friction coefficients for a 
single-mode cavity. Such an analysis serves as a ground 
to identify the physical processes underlying the cavity- 
cooling. In addition, it is interesting to compare the re- 
sults with those obtained in a similar system |7) with 
the cavity mode being driven instead of the atom. 

There were several attempts to interpret the cooling in 
terms of simple physical processes, such as two-photon 
Doppler effect Q, or Sisyphus-effect in the dressed state 
basis 0. In what follows, we propose to explore the 
friction mechanism by systematically varying the param- 
eters, through relatively simple limiting cases to a final, 
quite general parameter setting in the strong-coupling 
cavity QED regime. In Figs. || and || the friction coeffi- 
cient is shown as a function of the detunings in a contour 
plot style for various characteristic values k and g of the 
cavity. 

Let us start in the bad-cavity limit with a negligibly 
small coupling parameter (k = IO7, g — 7/2). This 
corresponds to the simple perturbative regime of cavity 
QED, where the role of the cavity is that it reshapes 
the radiative environment of the atom and increases the 
spontaneous emission rate at the cavity frequency. In 
our one-dimensional model we are interested in mechan- 
ical forces due to cavity photons. As these photons leak 
out of the cavity fast, such that no dynamics occur on a 
time scale longer than k , the atom-cavity field interac- 
tion can be basically considered as scattering. Although 
the atomic dipole is linear in the low-saturation regime, 
inelastic scattering can happen due to the CM motion 
that is able to compensate for the energy difference. For 
luc > ujp, the spontaneous emission being favoured at 
the cavity frequency higher than the incoming one, the 
scattering is accompanied by a loss of the kinetic energy, 
i.e. cooling. This is the origin of the cooling region for 
Ac < in Fig. ||a, and in turn, the same mechanism 
acts reversly leading to heating for Ac > 0. Note that 
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FIG. 2: Contour plot of the friction coefficient in the bad- 
cavity regime (k = 107) f° r various coupling constants, (a) 
g = 7/2, (b) g = 37/2, (c) g — 37, linking the closely in- 
dependent atom and cavity field system to the dressed atom 
one. Solid contour lines indicate cooling (/3 < 0), dashed ones 
heating (/3 > 0) regions. 



this simple energy conservation argument is also in the 
heart of the interpretation given in p] for the two-photon 
Doppler effect. 

Next, let us keep k large and increase gradually the 
coupling constant g. Figures ||b and c correspond to 
g = I.57 and g — 37, respectively. Two sharp peaks 
emerge at « 0. As k is still the far largest param- 
eter, the scattering picture applies. However, instead of 
the bare atom, the strongly coupled atom-cavity system 
has to be taken into account. As a consequence, rela- 
tively close to the pumping frequency lo p , the spectrum 
now exhibits the doublet of the first excited dressed-state 
manifold. Scattering of a pump photon into a cavity one 
may happen now via both intermediate states. In the 
case of Ac < and A^ « 0, for example, which is the 
lower half plane of the plots |p and c, the upper dressed 
state |+) is at about the cavity frequency ll>c with a width 
of k. The upper level therefore amounts to cavity photons 
of frequency about ojc, which corresponds to the process 
we described just before and is responsible for the broad 
cooling peak in the background. For large enough cou- 
pling constant, however, the lower state |— ) is contami- 
nated with a non-negligible amount of \g, I) component, 



its weight is proportional to g 2 /(ujc — oja) 2 ■ Pumping 
photons can then be transmitted into the cavity by ex- 
citing the lower dressed state at the frequency of about 
loa + 2g 2 /(uJA — t^c)- For A^ sa this state is more res- 
onant with the pump and the corresponding scattering 
channel becomes dominant. That is, one gets the heating 
peak for Ac < and A^ > 0, whose width is approxi- 
mately 7, the one of the lower dressed level. For A^ < 
the driving field is tuned below the lower dressed level 
and thus scattering via this state yields cooling again. 
Note the displacement of the heating peak with respect 
to the axis A^ = which is due to the cavity-vacuum 
induced lightshift, i. e., the lower dressed-state resonance 
and the atomic resonance do not coincide. 

Regardless of the value of the coupling constant, as 
long as <? <§; k, the field adiabatically follows the atomic 
dynamics, as is the case for the parameter settings chosen 
for Fig. [|. The damping effect cannot be attributed to a 
non-conservative, time-delayed field dynamics. This kind 
of behaviour occurs when decreasing the cavity linewidth 
k. Then, instead of being a passive element with specific 
mode density, it becomes necessary to include the cavity 
field as a dynamical component of the system. For longer 
spontaneous lifetimes, the system spends some time in 
the dressed states and the slow atom moves in a potential 
associated with the sinusoidally varying dressed levels. 

In Figs. ||a-c, the case of k — 57, 27, and 7 are shown, 
respectively. There appears a pair of peaks with hyper- 
bolic shape along AcA^ ss g 2 , which is cooling for nega- 
tive detunings Ac, A^ < 0, and heating for positive ones 
Ac, A a > 0. This cooling (heating) region corresponds 
to the situation where the lower (upper) dressed state 
is resonantly pumped at antinodes (maximum coupling) 
and thus the slowly moving atom has to, on average, 
climb up potential hills (descend to potential valleys). 

Finally, there are two peaks at at Ac ~ that can be 
attributed to the effect of the Doppler-shift in the corre- 
lated atom-field dynamics. Compared to the free-space 
Doppler cooling case, the preferential direction appears in 
the emission rather than in the absorption. For A^ > 
and Ac ~ 0, the atom gets closer to resonance with the 
copropagating component of the standing- wave field, i.e., 
following an emission the atom is more likely to impart 
a recoil in the direction opposite to its velocity. 

One can expect that in the low-saturation regime, 
where the Heisenberg-Langevin equations are linearized, 
it makes no difference which component of the coupled 
system is being pumped. Indeed, this intuition is re- 
flected in the mathematical form of the results. Wc 
checked that a systematic exchange of the parameters 
(A^, 7) <-> (Ac, k) in the expressions (|i~8l ) and (|i~f^ ) re- 
produces the results obtained for the cavity-driven case 
in Ref. @ . This symmetry reveals that the roles of the 
two oscillators, the field mode and the atomic dipole, are 
interchanged. Accordingly, the map of the friction coef- 
ficient in Fig. [|c is similar to the Fig. 3 of Ref. with 
the detuning axes exchanged (reflection to the diagonal) . 
In principle, for k = 7 there is a one-to-one correspon- 
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FIG. 3: Contour plot of the friction coefficient in the dressed- 
atom regime (g = 37) for decreasing cavity linewidth, (a) k = 
57, (b) k = 27, (c) k = 7, linking the adiabatic field dynamics 
regime to the one including time-delay, and the accompanying 
non-conservative cavity forces. 



dence between the two systems, i.e. a given dynamics 
of the atom-driving configuration can be mimicked, with 
exchanged detunings A c and A^, in the cavity-driving 
one. However, for a fixed setting of the detunings, the ac- 
companying cooling (or heating) mechanism is different, 
which leads to an essential modification of the relevant 
physical quantities. This deviation becomes of impor- 
tance when there is an additional constraint concerning 
the detunings, for example, has to be very large to 
keep the spontaneous photon scattering low for molecule 
cooling. The more detailed analysis of this comparison is 
relegated into the next section, where also other numeri- 
cal examples for the characteristic statistical features are 
presented. 



V. TEMPERATURE AND COOLING TIMES IN 
A MULTIMODE CAVITY 



In this section we will study thermodynamic properties 
of the system. Most importantly, we calculate the steady- 
state temperature which can be estimated by the ratio of 
the spatially averaged diffusion and friction coefficients. 
Localization effects were proven 13 to be important for 



cavity fields with higher intensity than the sub-photonic 
fields occuring in the present scenario. Hence, uniform 
position distribution of the atom can safely be assumed. 
In this approximation, the temperature becomes inde- 
pendent of the pumping strength rj. The other important 
feature is the time scale of reaching the given steady-state 
temperature. The so-called cooling time can be consid- 
ered to be the inverse of the friction coefficient (3. How- 
ever, it depends on the pumping strength 77 which is quite 
arbitrary, of course, within the limit of not to excite the 
atomic dipole too much. The interesting, pumping inde- 
pendent quantity, in fact, is the number of spontaneously 
scattered photons during the cooling time /3 _1 . Low 
number of spontaneous emission means efficient cooling, 
where the cooling time scale due to the cavity dissipation 
channel is short enough compared to the spontaneous 
scattering into lateral modes. In the following analysis 
we include the scaling of these two quantities on the num- 
ber of degenerate modes of the cavity. 



A. Effective mode approach 

We will study this effect first in a simplified geom- 
etry when the spatial variation of the different modes 
in the degenerate manifold is closely uniform. This can 
happen, for example, with a piece of coated waveguide 
where many quasi-degenerate sinusoidal modes can be 
found within the atomic spectral linewidth with slightly 
different wave numbers. In the following, we will con- 
sider another example, the first few transverse modes in 
a confocal cavity. The transverse modes with even in- 
dex are exactly degenerate. The corresponding mode 
functions, the Hermite-Gaussian modes, are known in 
the paraxial approximation. To be conform with it, the 
maximum transverse mode indices we can take into ac- 
count are limited by (n + m)A/Z cav <C 1, where / cav is 
the cavity length. This limitation also implies that the 
mode functions can be simplified around the cavity cen- 
ter z = 0. First, the Guoy phase term can be omit- 
ted, and the variation of the longitudinal wave number 
k2n, 2m = ^o,o — 2 (2b + 2m + l)// cav can be neglected, i. e. 
k2n,2m ~ ^o,o = k. Accordingly, the mode function along 
the cavity axis can be approximated by the simple cosine 
function cos (fez). Second, only the leading term of the 
derivative d z g(z), expanded into a power series of A/Z cav , 
must be kept, which is proportional to fcsin(Az). In a 
region close to the cavity axis the mode functions thus 
overlap and form an effective mode with enhanced cou- 
pling to the atom. It is an interesting problem to move 
out from this limit into a situation where the different 
modes have highly varying derivatives in space, that one 
has to study in a three dimensional context with the ex- 
act mode functions. 

In the simple example described above the presence of 
many degenerate modes can be incorporated in an effec- 
tive coupling constant g eH , a concept already used in Ref. 
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The enhancement factor is 



eff/ff 



(2iV + l)!! 
(27V)!! 



(21) 



where 2N is the maximum index taken into account, that 
is we consider a total number of modes M = (N + l) 2 in 
the dynamics. The effective coupling constant g Blt grows 
closely as a linear function of N, which indicates that or- 
ders of magnitude can be gained in the coupling strength. 
The unphysical divergence for large N stems from the ex- 
tension of the mode functions obtained within the parax- 
ial approximation to high indices of N. In practice, the 
effective g could be measured and then an effective num- 
ber of modes can be determined. 

Let us now see how the steady-state properties of the 
system depend on the number of modes in the confocal 
cavity example, i.e., on the effective coupling constant. 
The steady-state temperature is plotted in Fig. |]a for 
the hyperbolic cooling regime with Ac, A. a < 0, and 
A c « 0. When varying g cl{ , the detunings Ac and A^ 
are rescaled such that their product is fixed at <j£ (f , and 
their difference is constant. The first of these conditions 
ensures that the lower dressed state is pumped resonantly 
at an antinode (minimum energy). The latter one means 
that only the pumping frequency u> p is to be varied, both 
the atomic uja and the degenerate mode frequency cue are 
fixed. For the plot we set A^ — Ac = loc — — — 5O7. 
Since A^ ~ — 5O7, i.e., the driving field is tuned to 
be very far from resonance. The cavity properties are 
partly determined by the parameter k. The solid curve 
in the figure corresponds to k = 7, the dashed one to 
k = 7/IO. The other relevant cavity parameter g etl is 
considered a variable, however, it is useful to define the 
single-mode coupling constant g. It is then set to g = 37 
and g = 37/IO, respectively, as if the cavity length were 
changed by a factor of 10, keeping the same mirror trans- 
mittivity. Having defined the single-mode coupling g, a 
discrete series of effective coupling constants is obtained 
for increasing number of modes, which is indicated by 
the points on the curves. The rapid drop in temperature 
obtained when the number of modes included slightly in- 
creases is the main benefit we can expect from the use of 
a multimode cavity. In both cases, k = 7 and k = 7/IO 
sub-Doppler temperatures can be achieved with a small 
number of modes involved. Note also that the curves 
converge fast, indicating the existence of a well-defined 
temperature insensitive to the divergence of the effective 
coupling constant g ot{ with increasing number of modes. 

In addition to the induced atomic dipole moment, the 
driving field yields a small, but not completely negligible 
atomic excitation. Hence, spontaneous photon scattering 
occurs with a rate of 27 (cr^a). It is an important quan- 
tity how many photons are scattered in this way during 
the characteristic time of cooling which is (2/3) -1 in our 
case. The principal goal is to restrict this number be- 
low one which implies that the scheme can be extended 
for cooling particles without closed pumping cycle. The 
Fig. Hb shows that the number of spontaneously scat- 
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FIG. 4: Steady-state temperature (a), and the number of 
spontaneously scattered photons in a cooling time (b) versus 
the effective coupling constant g a tc. Dashed line corresponds 
to k = 7, the solid one to n — 7/IO. The crosses put on the 
curves indicate the coupling constant corresponding to maxi- 
mum transverse indices of 2N = 0, 2, 4, 8, with starting a 
single-mode coupling constant of g = 3y for the dashed curve 
of k = 7, and g = 37/IO for the solid curve of n — 7/IO. 
For this latter the single-mode coupling (maximum index 0) 
is missing from the plotted range of temperatures. 



tered photons in a cooling time, iV ph = j(a' f a}/f3, can 
decrease below the limit of one photon for large enough 
g eS , corresponding to about N = 64. 

Both quantities plotted in Fig. ^ are independent of the 
pumping strength. The cooling time itself depends on it. 
However, without speeding a driving field intensity, one 
can deduce numerical values of the cooling time from 
the Fig. |^b, provided the saturation is kept fixed. For 
example, at a saturation (cr^a) « 0.1, and for Rb with 
7 = 20//XS, the cooling time is N ph /4 in units of /xs. 

The result exhibited in Fig. ^ suggests that smaller k 
provides better performances in terms of temperature, 
cooling time. On the other hand, as it was pointed out 
in [j5], |Ig| |, the velocity capture range is limited by kv < k. 

Finally, let us return to the problem already addressed 
in the last section. Fig. || presents an additional curve 
(dotted line) that corresponds to the same parameter set- 
ting as the one belonging to the solid line (k = 7/IO), 
however, the single-mode cavity field is being pumped 



9 



instead of the atom. Whatever component is driven ex- 
ternally, the pumping field, by construction of the detun- 
ings, is resonant with the lower dressed state at an antin- 
ode, and this analogy makes the comparison justified. It 
is apparent that the temperature in the atom-driven case 
is lower. The curve belonging to the cavity-driven case 
for k — 7 does not even fit in the plotted range, the dif- 
ference with respect to the dashed line is much larger. 
Although a simple transformation connects the results of 
the atom- and cavity-driven cases, as this example illus- 
trates it, a significant physical difference can occur. The 
origin is that the detunings were the same, is large 
and Ac is small, which breaks the symmetry between 
the two systems based on the interchange of the above 
detunings. 



B. Beyond the effective mode approach 

The effective mode approach could be used in the pre- 
vious analysis because all the relevant modes closely over- 
lapped in the region of interest, i. e., around the cavity 
center. Accordingly, the system is reminiscent of a single- 
mode one with enlarged coupling constant g aB . By con- 
trast, when one moves away from the center, but still 
on the axis, the cosine-like mode functions with varying 
wavenumbers undergo a dephasing. One immediate con- 
sequence is that the friction force, proportional with the 
gradient of the mode function, does not vanish in any 
point. Furthermore, on inspecting the general solutions 
( |l8| ) and (|l9|), one can notice that the friction coefficient 
/? is proportional to the square of the sum d z J2 9ki z )i 
while the diffusion coefficient, for Ac = contains only 
the sum XX^Sfc) 2 - Although the determinant \D\ 2 ap- 
pears also on different powers in the denominator, it is 
clear that for certain parameter settings the two coeffi- 
cients can scale in a different way with the number of 
modes. This gives rise to the possibility to get an in- 
terferometric enhancement of the cooling by a collective 
effect of the modes and leads us to conjecture that the 
steady-state temperature can be lower in other positions 
than the cavity center. 

To check this expectation, we calculated the tempera- 
ture as a function of the position in the cavity by using 
the mode functions cos{(fcz — (2m + 2n + l)atan(z/zo)} 
with transverse indices n, m = 0, 1,N } where the Guoy- 
phase term is responsible for shifting the modes out of 
phase. The length scale of the dephasing can be esti- 
mated by z// cav ~ 7r/4(2A +1), i. e., at this distance 
from the center the Guoy-phase shift is about ir/2. For 
each position, (i) we perform again the averaging of the 
diffusion and friction coefficients over a wavelength, and 
(ii) we redefine the detunings A a, Ac such that the 
lower-dressed state corresponding to the local coupling 
constant be resonantly pumped at antinodes. The result 
is plotted in figure |^. The temperature initially drops 
as the atom moves out from the center. The higher in- 
dices are taken into account, the faster the initial drop 
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FIG. 5: Steady-state temperature as a function of the posi- 
tion in the cavity for k = 7/IO, g = 37/IO. The temperature 
in the center is used as a unit, therefore purely the competi- 
tive effects of the dephasing and the decrease of the coupling 
strength when moving away from the waist can be observed. 
The curves correspond 0, 2, 4, 8, 16 and 32 as maximum index 
taken into account. 



happens, which suggests that the underlying reason is in- 
deed the dephasing of the cosine modes. The estimated 
length scale shows a good agreement with the numeri- 
cally obtained results for various maximum indices 27V. 
The reduction of about 20% in the temperature can be 
attributed to a collective, interference-like effect of the 
multimode field. For large distances from the center, af- 
ter the dephasing is completed, the temperature grows 
slowly again, exhibiting the effect of the decreasing cou- 
pling constant. 

VI. CONCLUSIONS 

The mechanical effects of a high-Q cavity on the 
center-of-mass motion of a coherently-driven neutral 
atom have been investigated. We calculated the diffu- 
sion coefficient and the friction force from a quantum 
model adapted to the strong atom-field coupling regime, 
hence the validity does not depend on any specific re- 
lationship of the parameters. The model is analogous 
to the one described in Refs. || pj, however, we consid- 
ered the system with the atom being externally pumped 
instead of the cavity field mode. Surprisingly, this differ- 
ence leads to important new features in the diffusion and 
damping process. Lower temperatures can be achieved in 
the present system. We pointed out that certain cooling 
mechanisms can be realized only in the good-cavity limit, 
i.e. 5 > ^,7, where the dressed-atom dynamics including 
Rabi oscillations becomes dominant. The corresponding 
parameter regimes are especially suited to applications, 
since large atomic detunings (red or blue) are allowed 
here. Furthermore, we found that drastic improvement 
in terms of low temperature and small number of spon- 
taneous scattering can be obtanied in a degenerate mul- 
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timode resonator. As a highlight of this possible bene- 
fit, we showed that the number of spontaneously emitted 
photons from the atom during the cooling time can be 
reduced to below one, which demonstrates the possibility 
of cooling molecules optically below the Doppler limit. 
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